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摘 要 本 文 基于 OpenMP 并 行 语言 , 采用 拼接 式 的 多 块 格子 Boltzmann 方法 研究 了 接触 压力 、 粗 燃 度 、 间 际 介 质 导热 
系数 对 接触 热 阻 的 影响 规律 . 结果 表明 : 铝 之 间 的 接触 热 阻 随 着 接触 压力 的 增加 而 降低 ， 随 着 粗糙 度 的 增加 而 增加 ; 当 块 


体 气 凝 胶 之 间 的 间隙 填充 空 
触 热 阻 是 填充 空气 的 大 约 50 倍 . 
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气 时 , 气 凝 胶 之 间 的 接触 热 阻 随 压 力 的 增加 变化 不 大 ; 而 当 气 
随 压 力 的 增 大 而 降低 ; 间隙 介质 的 气体 导热 对 于 接触 热 阻 起 着 关键 的 作用 , 在 低 接触 压力 时 , 气 凝 胶 间 除 趋 于 真空 
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Abstract In the present paper, a patched-type multi-block lattice Boltzmann method is adopted to 


study the influences of contact pressure, roughness, gas thermal conductivity of interstice media on 


the thermal contact resistance. The results show that thermal contact resistance of aluminum pairs 


decreases with the increasing contact pressure, while increases with roughness. Besides, the thermal 


contact resistance of bulk aerogel pairs changes slightly with the increase of contact pressure when 


the interstices are filled with air, while significantly decreases if the interstice is approaching vacuum. 


In addition, the gas conduction of interstice media plays a key role on the thermal contact resistance. 


At low contact pressure, the thermal contact resistance of the vacuum gap is about 50 times of that 


of the air gap. 
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真实 表面 都 是 粗糙 表面 ， 因 而 当 两 材料 相互 接 
触 时 ， 只 有 很 少 接触 点 会 接触 上 ， 这 导致 了 接触 热 
阻 . 接触 热 阻 在 热 防 护 系 统 、 电 子 设备 中 扮演 着 重 
要 的 作用 。 从 1950 年 开始 ， 有 了 很 多 预测 接触 热 
阻 大 小 的 研究 工作 . 接触 热 阻 的 大 小 取决 于 材料 表 
面 形 貌 ， 接 触 变 形 机 制 、 及 传 热 分 析 。 对 于 表面 形 
貌 ，Greenwood 和 Williamson 出 发 现 , 粗糙 表面 的 
轮廓 高 度 近 似 服从 高 斯 分 布 . 在 这 个 模型 中 , 用 一 些 
统计 参数 , 如 粗糙 度 、 曲 率 , 来 描述 粗糙 表面 . 然而 ， 
这 些 统计 参数 和 测量 仪器 的 分 辨 紊 是 相关 的 ， 且 由 
于 粗 烟 表 面 形 貌 的 多 尺度 特性 ， 并 不 能 唯一 描述 表 
面 形 貌 . Majumdar 和 Tien P 采用 分 形 理论 来 描述 
粗 烟 表面 的 形 貌 ,在 这 之 后 ,基于 分 形 理论 的 W-M 


收 稿 日 期 2016-12-18; 修订 日 期 2017-02-20 


基金 项 目 : 国家 自然 科学 基金 国际 合作 重点 基金 (No.51320105004) 


multi-block; lattice Boltzmann method; thermal contact resistance; interstice media 


函数 日 被 广泛 用 于 生成 粗糙 表面 . 对 于 接触 变形 机 
制 ， 存 在 塑性 变形 、 弹 性 变形 ， 弹 塑性 变形 三 种 模 
型 . 弹性 变形 理论 基于 Hertz 理论 , 而 塑性 变形 认为 
接触 压力 超过 弹性 极限 时 发 生 的 都 是 塑性 变形 。 本 
文采 用 塑性 变形 理论 来 分 析 接 触 变形 机 制 . 

在 材料 的 接触 面 上 ， 存 在 三 种 传 热 途 径 : 1) 通 
过 接触 点 的 导热 ; 2) 通过 间 质 介质 的 导热 ; 3) 通过 
辐射 四。 在 以 往 的 研究 中 ,很 多 学 者 忽略 了 间 质 介 
质 (气体 ) 导热 和 辐射 的 贡献 553. 因而 本 文 的 目的 
是 为 了 考虑 气体 导热 和 辐射 对 接触 热 阻 的 影响 。 描 
述 具 有 一 定 分 辩 率 的 粗糙 表面 需要 大 量 的 网 格 ， 为 
节省 计算 资源 ,本 文采 用 如 图 1 所 示 多 块 网 格 的 方 
法 来 摘 述 接触 界面 . 在 接触 区 域 采 用 细 网 格 , 而 对 于 
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接触 材料 则 采用 粗 网 格 。 在 粗细 网 格 的 界面 上 需 保 
证 热流 、 温 度 等 物理 量 连 续 . 格子 Boltzmann 方法 
具有 易于 并 行 计算 的 特点 ， 因 而 本 文 基 于 OpenMP 
并 行 语言 ， 采 用 格子 Boltzmann 方法 用 于 数值 预测 
接触 热 阻 . 在 验证 了 多 块 Boltzmann 方法 程序 的 准 
确 性 后 , 研究 了 接触 压力 、 粗 糙 度 、 接 触 材 料 的 导热 
系数 、 间 际 介 质 导 热 系 数 、 温度 、 以 及 辐射 对 接触 热 
阻 的 影响 规律 . 


1 数值 方法 
计算 区 域 如 图 1 所 示 : 


Tn 


block 1 


coarse grid material 1 


block 2 
fine grid 


block 3 


coarse grid 


material 2 


7 
图 1 接触 热 阻 的 示意 图 


Fig. 1 Schematic diagram of thermal contact resistance 


对 于 整个 计算 区 域 , 能 基 传 输 方程 如 下 所 示 : 


oT 0/A oT 
a 下 
Ot Ti ( 闪 疾 )+5 WD 
边界 条 件 为 : 
oT 
了 | 一石 Tl,o=T2; 和 a- =0 (2) 
! l ” oz | -ao 


其 中 ,了 是 温度 , 入 是 导热 系数 , pcy 是 体积 比热容 , 5 
是 源 项 . 如 果 没有 内 热源 , 则 5=0. 
1.1 一 般 格 子 Boltzmann 方法 
在 格子 Boltzmann 方法 ， 对 应 式 (1) 的 温度 分 
布 演化 方程 如 下 四 
fa(zt + eadt,t+6t) — falz,t) = 


1 (3) 
=(falz,t) — fa (7,t)) + dtwaS 


其 中 , z 代表 位 置 , t 是 时 间 , 5t 是 时 间 间 隔 , 大 是 
温度 分 布 函数 ， 大 4 是 平衡 态 温度 分 布 函 数 , 定义 为 


(D3Q7 模型 ) 门 ; 
f=waT, a=0—6 (4) 


其 中 ， Wa =1/7., Ea 一 是 离散 速度 : 


01 -1 0 0 0 0 
ea=|100 0 1 -1 0 0 |c (5) 
0 0 8 0 1 =<1 
热 扩 散 系 数 和 松弛 时 间 的 关系 : 
2 1 
a= 3 (3) (6) 
当地 温度 可 由 下 式 计算 得 到 : 
?= 0 
而 当地 热流 密度 四 ， 
gq 二 —pcp— = > eazfa (8) 


1.2 多 块 格子 Boltzmann 方法 
图 2 所 示 的 是 两 块 网 格 系统 的 示意 图 ， 用 于 说 
明 界面 上 信息 的 传递 过 程 . 网 格 比例 : 


m= 6xc/07T8 = te/otr (9) 


coarse grid 


fine grid 


图 2 两 块 网 格 系统 
Fig. 2 Two-block grid system 


可 以 推导 出 ， 网 格 上 的 信息 量 从 细 网 格 传递 到 
粗 网 格 的 过 程 可 由 下 式 获得 和 


+ peq ,Tl A + eq 
2 tp fy) (10) 


其 中 , f++ 是 碰撞 后 的 温度 分 布 函数 , 下 标 c 代表 的 
是 粗 网 格 ,而 下 标 了 代表 的 是 细 网 格 . 从 粗 网 格 伟 
递 到 细 网 格 的 过 程 可 由 下 式 计算 得 到 ， 

ts = +m 


TH lM /pt peq 
二 入 (fas 54) (11) 
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式 (10) 和 式 (11) 可 以 确保 热流 在 粗细 网 格 界面 
上 连续 . 

值得 注意 的 是 , 在 界面 上 黑色 实体 格 点 。 的 信 
息 需 要 插值 得 到 ， 我 们 采用 三 次 样 条 插值 获得 . 粗 
细 网 格 的 步 进 时 间 步 长 也 不 一 样 ， 时 间 上 的 插值 采 
用 三 点 拉 格 朗 日 来 获得 . 
1.3 接触 热 阻 

接触 热 阻 模型 的 示意 图 如 图 1 所 示 。 粗 网 格 
100x100x100, 细 网 格 800x800x30。 为 保证 热流 在 
气 固 界面 上 热流 连续 , 还 应 使 : 

(pep)p = (pcp)s (12) 

这 样 的 假设 , 当 达 到 稳定 状态 时 ,并 不 影响 温度 场 . 
对 于 图 1 所 示 的 计算 区 域 ， 上 下 给 定 等 壁 温 ， 边 界 


绝热 . 则 当 温度 场 稳定 时 ， 接 触 热 阻 可 由 下 式 计 算 


4 租 到 [10], 
得 到 A 


其 中 , At 是 接触 界面 的 温差 , q 是 热流 密度 . 粗糙 表 
面 采用 W-M 函数 获得 : 


z(z,y)=L (中 ( 旦 ) 1/2 


之， YO { cosgmn— 


L 
oon (tm () -加 ) +6mn| 
接触 压力 由 下 式 计算 得 到 办 
全 元 (15) 


(a) p=0.116 MPa 


(b) p=5.179 MPa 
图 3 接触 表面 对 


Fig. 3 Contact surface pairs 


生成 的 粗糙 表面 对 如 图 3 所 示 。 
2 模型 验证 

图 4(a) 给 出 了 接触 热 阻 的 简单 计算 模型 ， 用 两 
个 固体 之 间 夹 着 一 层 空 气 来 模拟 接触 热 阻 . 间隙 距 
离 为 5.2 um, 间 际 介质 导热 系数 0.026 W/(m'K), 固 
体 导 热 系数 1 W/(m:K). 粗 网 格 11x11x11, 细 网 格 
21x21x21. 用 多 块 格子 Boltzmann 方法 预测 得 到 的 
接触 热 阻 为 2x10 忆 m2.K/W, 和 理论 预测 值 一 样 说 
明了 该 方法 的 准确 性 . 从 图 4(b) 可 以 看 出 ,热流 密 
度 在 粗细 网 格 界面 上 连续 . 
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不 多 ， 因 而 当 间 际 介 质 填充 空气 时 ， 气 凝 胶 固体 之 
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(b) 沿 着 z 轴 的 热流 密度 
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图 4 接触 热 阻 的 简单 模型 


Fig. 4 _ Simple model of thermal contact resistance 


3 结果 与 讨论 
3.1 接触 压力 的 影响 

图 5 给 出 了 铝 接触 热 阻 随 着 接触 压力 的 变化 关 
系 . 铝 的 平均 表面 粗糙 度 是 R。=1.95 km。 从 图 中 可 
以 看 出 ， 随 着 接触 压力 的 增加 ， 接 触 热 阻 先是 急剧 
降低 而 后 下 降 得 平 级， 接触 压力 的 增加 导致 了 接触 
间 陈 的 减 小 ， 而 间 阶 介质 的 导热 系数 比 接触 材料 的 
导热 系数 低 。 因 而 随 着 接触 压力 的 增加 ， 材 料 之 问 
的 接触 热 阻 会 降低 . 


TCR/m:K:W! 
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图 5 接触 热 阻 随 压力 的 变化 


Fig. 5 Pressure dependence of thermal contact resistance 


3.2 间隙 介质 导热 系数 的 影响 

气 凝 胶 常 用 在 飞行 器 的 热 防 护 系统 中 . 本 文 研 
究 了 块 材 气 凝 胶 之 间 的 接触 热 阻 . 气 凝 的 导热 系数 
设 为 0.03 W/(mnK)。 和 气 凝 胶 表 面 的 平均 粗糙 度 为 
有 Rs=1.95 um。 图 6 中 给 出 了 间隙 介质 导热 系数 对 
气 凝 胶 接触 热 阻 的 影响 规律 . 由 于 气 凝 胶 固 体 本 身 
的 导热 系数 和 空气 的 导热 系数 (0.026 W/(m'K)) 差 
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间 的 接触 热 阻 随 压力 增 大 的 下 降 趋 势 并 不 明显 . 而 
当 间 隙 介质 导热 系数 2.6x10“W/(mnK) 时 ， 气 凝 
胶 固体 之 间 的 接触 热 阻 随 着 压力 的 增 大 而 增 大 . 因 
而 , 间 院 介质 的 导热 系数 的 对 接触 热 阻 的 影响 很 大 ，. 
在 低 接 触 压力 时 , 间隙 介质 的 导热 系数 为 2.6x10 
W/(m-K) 的 接触 热 阻 是 间隙 介质 导热 系数 2.6x10 
W/(m-K) 时 的 大 约 50 倍 . 


—e— aerogel, 4,=2.6x10 ? W/(m-K) 


~® acrogel, 4,=2.6x10" W/(m-K) 
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图 6 间隙 介质 导热 系数 对 接触 热 阻 的 影响 


Fig. 6 Gap medium influence on thermal contact resistance 


3.3 表面 粗糙 度 的 影响 

表面 形 貌 会 对 接触 热 阻 产生 很 大 的 影响 . 这 里 ， 
我 们 研究 了 铝 材料 表面 的 粗糙 度 对 接触 热 阻 的 影响 。 
图 7 给 出 了 表面 粗糙 度 对 接触 热 阻 的 影响 规律 . 这 
里 的 表面 粗糙 度 是 平方 根 平均 粗糙 度 . 从 图 中 可 以 
看 出 ， 接 触 热 阻 随 着 粗糙 度 的 增加 而 增加 。 这 是 因 
为 随 着 粗糙 度 的 增加 ， 上 下 接触 表面 之 间 的 平均 距 
离 增加 , 而 间隙 介质 之 间 的 导热 热 阻 是 主要 热 阻 . 因 
而 随 着 表面 粗糙 度 的 增加 ,接触 热 阻 也 会 增加 . 


1033 


TCR/m2. 天 .W-! 


R,/um 
图 7 表面 粗糙 度 对 接触 热 阻 的 影响 


Fig. 7 Influence of roughness on thermal contact resistance 
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胶 之 间 的 接触 热 阻 随 压力 的 增加 变化 不 大 ; 而 当 气 [6] Mishra S C, Poonia H, Vernekar RR, et al. Lattice Boltz- 
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作用 ， 在 低 接触 压力 时 ， 气 凝 胶 间 际 趋 于 真空 时 的 [7] Xuan Y M, Zhao K, Li Q. Investigation on Mass Diffu- 
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